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1. Introduction 

The perturbativc expansion of the off-shell amplitudes for a supersymmetric gauge the- 
ory may be plagued by infrared divergences for a number of different motivations. One 
possible source of infinities is given by the presence of positive mass-dimension couplings 
associated to massless fields in superrenormalizable theories. Canonical Yang-Mills the- 
ory coupled to massless matter in three dimensions turns out to be a good playground 
to study these phenomena [1] . In this case, going high in the order of the dimensionful 
coupling in the expansion of a given amplitude, for dimensional reasons one obtains 
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high powers of external momenta in the denominator. If the fields are massless, upon 
inserting those amplitudes in higher order graphs, IR divergences inevitably show up 
in the Euclidean integrals. 

Considering theories with only dimensionless couplings greatly improves the sit- 
uation. A direct power counting argument shows that in four-dimensions classical 
marginality is a sufficient condition to exclude the presence of IR infinities in ampli- 
tudes with generic external momenta [2]. Even if this argument may be generalized 
to marginal three dimensional theories, it only applies to component field formulations 
and cannot be applied to superspace formulations as we will show throughout this work. 

In fact, a potential source of IR divergences has to be considered as soon as the 
computations are performed using supergraph techniques in supersymmetric gauge the- 
ories. While proving to be an efficient method to compute perturbative corrections, 
superspace algebra comes with additional infrared issues due to the peculiar nature of 
the gauge supcrficld propagator. As it is clearly described in [3] in the case of four- 
dimensional Super- Yang- Mills theories, the appearance of infrared infinities can be as- 
cribed to the presence of (corrected) vector lines in loop diagrams. With a canonically 
gauge-fixed action and omitting the color structure, the gauge superfield propagator in 
momentum space can be written as 

l±S^^^5\e-e') (1.1) 

where a is the gauge- fixing parameter and Vq — —^{D'^D^ + D'^D'^) is the super- 
spin zero projector. Recall that this operator satisfies Vq + V1/2 = l, where V1/2 — 
^D"D^Da is the superspin 1/2 projector. It is clear that, already at the one-loop 
level, the Fermi-Feynman gauge a = 1 is the only infrared safe choice. On the other 
hand radiative corrections, being governed by Slavnov- Taylor identities, come with the 
transverse structure 7^1/2, thus reintroducing the infrared dangerous part in the propa- 
gator. Therefore, unless a way is found to perturbatively maintain the Fermi-Feynman 
form of the propagator, IR divergences will show up again starting from two-loop order. 

An explicit prescription to cure the IR divergences in the case of four dimensional 
Super- Yang-Mills theory has been given in [4]. The main idea is to introduce a non- 
local gauge fixing term and renormalize the gauge fixing parameter to preserve the tree 
level structure of the Fermi-Feynman propagator. The prescription presented in [4] is 
strongly based on the nature of the model (gauge sector of Yang-Mills type) and on 
the space-time dimension. 

The aim of this paper is to study the infrared behaviour of supergraph amplitudes in 
the case of marginal Chern-Simons- matter systems in three dimensions described in the 
J\f — 2 superspace formalism. These models can be treated in strong analogy with 4d 
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Yang- Mills theories while exhibiting a completely different gauge structure. Due to the 
recent interest in CS- matter models in the context of AdS/CFT we start our analysis 
in the specific case of ABJM theory [5]. By directly computing Green's functions up 
to two-loop order, we show that infrared divergences appear in the amplitudes but can 
be seen as a gauge artefact of the formahsm. At first, we gauge fix the ABJM action in 
a canonical way (a-gauge) and show that, in analogy with the four-dimensional SYM 
case, IR infinities show up when a corrected gauge vector propagator is inserted in loop 
amplitudes. By direct inspection of the dependence of the infrared singularities on the 
gauge fixing parameter «, we conclude that there is no suitable choice for the latter 
that both eliminate the divergences and preserve hermiticity of the action. 

To solve this problem, we slightly revise the prescription of [4] introducing a set of 
non-canonical gauge fixing terms (77-gauge) . This new set of gauges has the virtue that 
it can be used, by perturbatively fine tuning the parameter r), to complete loop by loop 
the transverse structure of the gauge vector propagator with the longitudinal part, thus 
improving its behaviour in the infrared. We will show how the infinities are canceled 
in this way by direct perturbative computations. The 77-gauge can hence be considered 
as a tool to consistently study the perturbative expansion of the amplitudes of the 
model without the presence of IR divergences. It's important to stress that infrared 
divergences, being an artefact of the superspace formalism, do not manifest themself in 
physical gauge invariant quantities. In this case the a- and 77-gauges produce coincident 
results. 

As a byproduct of our analysis, we explicitly compute the finite expression in a 
general gauge for the two-loop propagator of the chiral superfield in ABJM. Moreover, 
we study a special vanishing external momenta limit of the two-loop vertex function 
of ABJM theory showing that it produces a finite result. Finally we comment on the 
extension of our results to a general perturbative order and to any classically marginal 
Chern-Simons-matter system. 

2. ABJM action and gauge- fixing 

To address the problem of IR divergences in three-dimensional Chern-Simons theories 
we restrict ourself to the specific case of ABJM model. This theory possesses remarkable 
properties such as extended supersymmetry and exact conformal invariance which will 
simplify the analysis of the infrared behaviour. We extend our results to more general 
CS theories in Section 4. 

At first, we set up our notations and quantize the theory introducing the standard 
gauge fixing procedure and an alternative one that will ensure the cancelation of the 
IR divergences in loop amplitudes. We will use the Af — 2 superspace formulation first 
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presented in [6] adapted to the notations of [7] (see Appendix A for further details). 
In Eudidean space, we quantize the theory with a path integral measure of the form 
jD(f) e^t"^!. ABJM theory can then be written a,s S — Scs + ^mat + Spot, where 



cs 



k 

An 



(fxire / dt tr 



VD' 



-.-tv 



tv 



(2.1) 



■'mat 



d'^xdW tr 



2TTi 



BD 



V 



2TTi 



tr(B^^ij-B^^s) + — 



(2.2) 
(2.3) 



The chiral superfields B^ and Aa (where a,b,c,d=i,2) transform in the (2, 1) and (1, 2) 
of the global SU{2) x SU (2) respectively. Moreover, they transform in the (N, N) and 
(N, N) of the gauge group Uk{N) x U-k{N), such that if explicit gauge group labeling 
is needed, the chiral superfields are B^, Al, with a, b,d,b = 1, .., A^. The gauge 

vector superfields V and V are in the adjoint representation of the groups and may be 
written either as V — T^Vi with i = 1, ■ ■ ■ , A^^ or with matrix labehng V^, VP. 

To quantize the theory we re-scale the vector gauge fields V \f^^ we 
choose in each gauge sector the gauge fixing functions F — D^V, F — D^V. The 
standard procedure in d = 3 is to introduce in the functional integral the factor: 



/ 



VfVfA{V)A'\V)exp 
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2^ 



d^xd'^etr{f f) exp 
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2Ei 



d'xd'etr (//)), 

(2.4) 



where 



A(V) = / 



d\ d~k S{F{V, kK)-f)6{F{V,KK)- /), 



with A the chiral superfield of the gauge transformation — )■ 



(2.5) 

and a a 



dimensionlcss parameter. Notice that, following [10], we are introducing a gauge av- 
eraging given by gaussian weights with chiral integrals of the form ~ ^^e^ . This 
is to be contrasted with the standard d — A, J\f — 1 superspace procedure where one 
averages with a non-chiral (whole superspace) gaussian weight of the form e^^^. The 
average produces the canonical quadratic gauge fixed action [8] 



c(«) 
^9f 



I / d^xd'^9 tr 
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(2.6) 
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such that, after inverting the kinetic operator we obtain the gauge field propagators in 
momentum space 



y4 _ 1 

d — ~ 

P 

with = — 0'). From now on we shall call this gauge fixing procedure as the 

"a-gauge". This propagator simphfies greatly when a — )■ (Landau Gauge). In the 
next Section we will see that, if we want to preserve the hermiticity of the gauge fixed 
action considering a and a as complex conjugates, then the infrared divergences cannot 
be canceled by a simple fine tuning of the gauge fixing parameter. 

To solve this problem, in analogy with [4], we propose a different gauge averaging 
procedure. We choose the same gauge fixing functions as before, but this time we 
introduce the following term in the functional integral: 



detM j VfVfA{V)A-\V)exp (^j d^xd^Oti {jMf 



(2i 



In this way we allow for a non-trivial gauge averaging by the insertion of the operator 
jCi. It is important to stress that as long as M. is field independent, the detTW factor 
appearing in the functional integral is irrelevant and there is no need to introduce 
Nielsen-Kallosh ghosts [9] in the action. In c? = 4 the A4. operator is dimensionless 
and one can simply choose A4. = constant. In three dimensions it has dimensions of 
length so that we may choose either a dimensionful constant or a non-local gauge fixing 
term. Our choice of this operator in momentum space is Xi{p) = ^ , where 77^^^ is 

a dimensionless function that contains e powers of p ^. The early introduction of the 
e = I — I regulator parameter has to be understood formally in the sense of dimensional 
reduction, that is, we will still perform D-algebra calculations in three dimensions and 
only at the end we will regularize Feynman integrals. More specifically, we will define 
r] as an odd power series in the 't Hoof coupling A = ^ with coefficients that we will 
conveniently choose. By choosing the same Gaussian measure (2.8) for both gauge 
sectors we obtain the gauge fixed action in momentum space 



SI" = - / j^d'^e tr 
f-^ 2 J (27r)3 



V{-p) I D^D, - ^Vo 1 V{p) 



tr 



V{-p) { D^D, + ^Vo ) V{p) 



(2.9) 



^This choice would introduce Nielsen-Kallosh ghosts in the action if the computations were per- 
formed using the background field method as in [4]. 
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Inverting the operators of the quadratic part of the gauge fixed action we obtain the 
gauge field propagators 



We shall call this gauge fixing procedure as the "?7-gauge" . We will show that allowing 
7] to be corrected order by order in the 't Hooft coupling ^, we will be able to cancel 
the infrared divergent parts in the amplitudes. To complete the gauge fixing procedure, 
we rewrite the A~^{y) factor in the path integral using the usual Fadeev-Popov h-c 
and h-c system of Grassmanian chiral superfield ghosts. In both of our gauge fixing 
choices the ghost action reads 



/" r- 1 All 

Sfp= / d^xd'^e tr^bc + cb+ -J — {b + b)[V,c + c 



-lc-fb-^^^ib + b)[V,c + B] + . (2.11) 

In Appendix B we detail some of the relevant Feynman rules for ABJM theory. 



3. Infrared behaviour of amplitudes in ABJM theory 

We would like now to understand the origin of the IR divergences in the perturbative 
expansion of the off-shell amplitudes. In order to do so, we compute the one-loop 
correction to the vector gauge superfield V in the a- and r^-gauges. Performing a 
direct two-loop computation of the (finite) corrections to the self-energy of the matter 
superfield and to the superpotential, we will see that IR infinities only arise when the 
one-loop corrected gauge propagator is inserted in loop diagrams. A suitable choice of 
the rj gauge fixing parameter will then cancel the divergences. Our explicit examples 
will be completed with an all loop analysis in Section 4. 

3.1 One- loop vector propagator 

The one loop corrected gauge vector field receives contributions from matter, ghost and 
gauge vector fields as we show in Figure 1. In the a-gauge these evaluate 
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Figure 1: One- loop self-energy corrections of the gauge vector superfields. 



Alile. = -f^'Aj V^^<-A:)Pv2n(A:) • (3.1) 



Here we are displaying the corrections to the V — V propagator as well as to the mixed 
V — V propagator (last line). The latter receives contributions only from the matter 
diagram of Figure 1. The corrections to the V — V propagator can be easily read from 
the V — V case. The definition of the functions can be found in Appendix C. 

Notice that the mixed V — V contribution is subleading in N. Consistently with the 
Slavnov- Taylor identities, the complete one-loop correction to — contains only the 
spin 1/2 projection of the gauge field and is given by 

^^v^ = ^ ((3 + ««)5X'iV + (1 - aa)«) J '^^^Gl^^.^ik^-^ Vl{-k)V,„V^,{k). 

(3.2) 

With the particular choice a = (Landau gauge) we exactly reproduce the results of 
[10] [11]. We notice that, looking at the leading part of the correction, the effect of 
working in a general a-gauge simply results in a positive constant shift Therefore 
we conclude that, if we want to preserve hermiticity of the action, there is no way to 
fully eliminate the one-loop correction by fine tuning the gauge-fixing parameter. 

Now we work out the correction in the //-gauge restricting the analysis to leading^ 
order in N . As anticipated, the idea is to make a perturbative expansion of the gauge- 
fixing parameter in the 't Hooft coupling 

Vlk)= S^fc)A + 0(A3), (3.3) 



^The following ideas may be also worked out at subleading order in N but we would have to add 
a mixed V — V gauge fixing term. 
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such that to order A we obtain the total correction of the propagator: 

^ (pyT72 (•^•4) 

It's easy to see that if we choose ^V(k) — i)ik'^)~'' exactly complete the 

transverse structure V1/2 with the longitudinal part Vq to obtain 



-671^11) A ^^2)7/2^.- (3.5) 



In the next Section we compute two-loop Green's functions and show that the improved 
IR behaviour of the //-gauge propagator in (3.5) is enough to cure the problem of IR 
infinities. 

3.2 Matter self-energy at two-loop order 
3.2.1 Landau gauge 

Working in the Landau gauge simplifies greatly the calculation since many diagrams 
can be discarded due to the form of the gauge vector propagator. It is easy to see 
that all one-loop corrections vanish with standard gauge averaging. In Figure 2 we 
display all non-vanishing self energy two loop quantum corrections of matter fields in 
this gauge. The blob in diagrams (c) and (d) represents the insertion of the full one 
loop correction to the gauge propagator. Any other potentially contributing diagram 
is zero either by D-algebra or by color symmetry. 




d. 

Figure 2: Two loop self-energy quantum corrections. 

Let us for example calculate with detail diagram b of Figure 2. Taking into account 
the possibility of having V-V, V-V and V-V internal lines and using color vertex factors 
as (B.7), it evaluates to 

n, = -1 (iv^ - 1) ' I d'e d'e' j^.tr {Ba{-p)b\p)) v,{9, &) (3.6) 
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with 

the D-algebra factor of the supergraph. As we mentioned before, we perform all D- 

algcbra manipulations in three dimensions and wc calculate the final Feynman integral 
in d dimensions. After the usual integration by parts we obtain an ultraviolet divergent 
contribution 

= {N' - 1) (f ) ' / ^ tr {B4-P)B^{P)) ^1,1)^1,1/2+.) iPr^'- (3.8) 

To obtain the contribution a from Figure 2 we need vertex factors (B.5) and (B.6). 
We get the UV divergent contribution Ha = 2Ili). 

Using the corrected vector propagator we obtain for graph d of Figure 2 an UV/IR 
divergent tadpole 

-3(^^-1) (^)Y;|^^t.(B,(_rt6^(rt) / (gj^. (3.9) 

and for graph c an infrared divergent contribution 

(3.10) 

Summing up a, b, c and d we obtain the cancelation of all UV divergent contributions 
and we are left with an IR divergent piece 

n„+,+e+. = 3 (AT^ - 1) ' I tr {B^{-p)B\p)) g,ip), (3.11) 

with Qdip) as given in the appendix. Finally, diagram e produces a finite correction 

n, = -2(Ar^-l) J^^^m-P)B^{P)) h (3.12) 

where If. is the factor obtained after closing the D- Algebra: 

d^k dH {k + pf{l + pf - kH^ + p\k + l+pf _ 1 



(27r)3(27r)3 k^k + p^ {k + I + p)-" P {I + p^ 64' ^^"^^^ 

To conclude, the sum of all contributions gives a finite and an infrared divergent piece 

Working in the Landau gauge, we explicitly see that infrared infinities are only given by 
graphs c and d, which correspond to insertion of the 1-loop corrected vector propagator. 
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3.2.2 a-gauge 

We now take the more general case a 7^ 0. Once again there are no one- loop matter 
corrections. The list of two-loop self energy contributions gets larger. Apart from those 
already displayed in Figure 2, which are modified by the more general a-dependent 
propagator, we have some additional contributions displayed in Figure 3. The new 




/• 9- h. 

Figure 3: Additional quantum corrections in the a-gauge. 



contributions produce additional UV, IR divergences and finite pieces. The sum of the 
original diagrams we had, with the modified a-dependent propagator gives 

X (-4aaq,_„G;,,i/2^,, (pY" + (3 + aa)&(p) - ^(1 + aa)'j . (3.15) 
And the contributions from the additional diagrams of Figure 3 

nj,,,, = {N'' - 1) y ^ {Ba(-p)B^(p)) 

X (4aa (P'r" + 3^aa) ■ (3.16) 

By summing up all the contributions, we find as expected that all UV a-dependent 
divergences cancel out. The finite piece we had already encountered in the Landau 
gauge is not modified, and the IR divergent piece gets shifted: 

n" = (iv^-i)(^)7^t^(^"^(-^)^^(^)) ((3 + - ^) • (3.17) 

From this we may conclude that if we only allow a hermitian gauge-fixed action, such 
that a is literally the complex conjugate of a, then it is not possible to choose a value 
of the gauge fixing parameter a such that the self-energy corrections are infra-red safe. 
This is a direct consequence of the fact that IR divergences are eventually produced only 
by corrected vector propagators. It's also important to notice that the finite correction 
to the propagator turns out to be gauge independent, even if the propagator itself is 
not a physical quantity. 
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3.2.3 //-gauge 

In the r^-gauge the vector superfield propagator is written as: 

m-p) V,%p)) = + ^Po) S^e,e'/X- (3.18) 

where 77^^-) is expanded as in (3.3). The first piece of the propagator gives rise to matter 
self energy diagrams starting from two loops with the same contributions as in the 
Landau gauge (see fig. 2) such that, for 1, it gives a finite and an infrared 

divergent piece of order (^)^ given by 

n^+ -+e = \-^) J tr {Ba{-p)b\p)) l^3g,{p) - -J . (3.19) 

The second part of the propagator produces one-loop corrections to matter self energy 
such that, if the gauge parameter is of order ^, the contribution is of order {^Y- In 
this case, two loop and higher corrections will contribute beyond so we do not 
consider them. The one-loop contributions are displayed in Figure 4. 




Figure 4: One loop corrections in the ?7-gauge. The small black squares in the gauge vector 
propagators should be intended as the rj dependent piece of the propagator. 

After straightforward D-algebra, and with the choice r/^^^ = — GttG^^ ^^{k'^)^'^X + O(A^) 
we obtain 

n?.,oop = - (^)7 {Ba{-p)b\p)) sg.ip). (3.20) 

From this we see that with the choice of r/^^-j that produces the IR improved gauge 
propagator, we cancel the infrared divergent part obtaining only the universal finite 
piece already computed in the a-gauge: 

We therefore conclude that the improved IR behaviour of the gauge propagator is 
sufficient to eliminate the presence of the unwanted divergences. In the next Section 
we further check this assertion computing at two-loop order the matter four-point 
Green's function. 
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3.3 Superpotential vertex corrections 

The set of all two-loop graphs which contribute to superpotential corrections to lead- 
ing order in N in the Landau gauge are depicted in Figure 5; any other potentially 
contributing 2-loop graph is zero due to color symmetry, supersymmetry or particular 
symmetries of the Feynman integrals involved. Notice that, since in the Landau gauge 
the one loop correction to the vertex is exactly zero, we can discard many diagrams at 
2-loops that contain the 1-loop diagram as a subdiagram. 




a. b. c. 




d. e. /■ 

Figure 5: All two loop quantum corrections of the superpotential. 

To simplify notation, whenever we put a V{- ■ ■) in front of the graph inside an 
equation, we mean the scalar graph with all the momenta in the numerator gener- 
ated after closing the D-algebra (we put on equal foot the V and V lines and only in 
color/flavor vertex factors will we consider the sign difference between their propaga- 
tors and couplings to matter). Else, in the absence of T> in front of the graph, we just 
mean the corresponding scalar Feynman integral. 

To leading order in N, all the two loop contributions produce a term proportional 
to the classical superpotential (no double traces are generated) given by 

r.l^. B] = (1^) ' C. I . . . (2.).'^ (p, + + p3 + p.) 

977-7 

— e^ce^^tr {B''{p,)Ad{P2)B'' {pMb{pa)) A(pi, ■■■ ,Pa) i = a,--- J, (3.22) 

where Cj is the vertex factor of graph i. 2^j(pi, ■ ■ ■ ,^4) is the Feynman integral which 
results after performing the D-algebra so as to eliminate all the d^6 integrals except 
for the last one which is used to transform the D-operators applied on the fields into 
external momenta by using that / = J d?9D'^{- ■ ■). The vertex factors for all graphs 
are Ca = ^, Cb = |, Cc = —3, Cd = 1, Cg = —1, Cj = 2. We will always consider 
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Pi as the 'north- western' momentum of the graph and name the consecutive momenta 
counter-clockwise. 

Let us start the computation of graph c of Figure 5 which is the one we expect to 
give an IR divergence. A not so straightforward calculation of the D-algebra gives 






(27r)^ (A;2)3/2+e(A; + p4)2(A;-p3)2 

\qM - l^diP^), (3.23) 

where we have written the Feynman integral in terms of a finite scalar integral and 
infrared divergent contributions. Once again we obtain infrared divergences when we 
attach a one-loop corrected gauge vector inside a loop. We expect that IR divergences 
in the superpotential are canceled by the exact same choice for rj we found to improve 
the gauge vector propagator infrared behaviour. This is in fact true: the one-loop graph 
with the ?7-dependent part of the gauge vector propagator gives 

-Q7rX{gd{ps) + Gd{P4)). (3.24) 

With the value of the gauge parameter we made before r] = —GnCf^^ i)(p^) " A + O [A^], 
this insertion produces a superpotential correction with the same structure as in (3.22). 
In this way, the sum of this graph with graph d which was also IR divergent gives 

-3(47rA)^p|^^<^j+47rAI?(^^^j =-^(47rA)2(p3+P4)^^, (3.25) 

which is finite. These are the only dangerous IR graphs contributing to the superpo- 
tential; the graphs which remain to be analyzed are all finite. To show this, we list the 
integrals resulting from D-algebra computations. 

The simplest graph is g: it has three possible channels of which two contribute to 
leading order in N. This factor of 2 is already taken into account in the vertex factor 
Cg. The D-algebra of this graph is simply 

d^k dH -{Pi+P2? ]_ ^3 2g) 




(27r)3 (27r)3 {k + pi + P2Y P (l - P3 - Pif 64' 



In all other diagrams, a sum over different distributions of internal lines has to be taken 
into account such that diagrams 6, c and e appear four times with different momentum 
distribution, while diagrams a, d appear eight times. 
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A straightforward calculation shows that 





2(P3+P4)^ M • (3-27) 



As mentioned before, this scalar integral is finite in three dimensions. For graph a we 
obtain the finite result 

' ' ' {2t^y{2t^y [k + p^y {k - p^y {k + ly {I - p^)^ p • i • j 

Notice that the presence of a three lined vertex is potentially dangerous, but the mo- 
menta in the numerator of the Feynman integral that we obtain through D-algebra 
guarantees finiteness. The same is true for graph d 

d^k dH -Tr (7^7.7^7,7^7^) {k + p,Y F jk + ly jk - PsYpIpJ 

F {k + piy {k - p^y {k + ly (/ + p^y p 

(3.29) 



(3.30) 




d^k dH -Tr (7^7,7^7^) k^ ^ 

{2TTy {2TTy p _ (k + i + p^y (/ - p^y p ' 

which are once again finite in three dimensions. 

3.3.1 A particular exceptional momenta configuration 

By using the r) gauge fixing, we showed in the last Section that it was possible to obtain 
an infrared safe function of the external momenta for the superpotential corrections. 
Moreover, it is clear that the sum of IPI graphs plus four-legged graphs with corrected 
legs (using self-energy corrections we derived before), is a physical gauge invariant 
quantity. Having found an universal finite value for the matter propagator correction 
we conclude that also the correction to the superpotential is (at least at two loops) 
gauge independent. We would like now to compute it for a special external momenta 
configuration ^. 

The calculation we are going to present here should be interpreted along the lines 
of [19]. In these papers, by means of direct computation of specific diagrams in four- 
dimensional supersymmetric models, it was shown that finite contributions may survive 
the limit of vanishing external momenta for the IPI vertex function as soon as massless 



^See [17] for the calculation of the effective action on a vector superficld background. 
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particles were present. Moreover, these contributions could break holomorphy in the 
coupling constants or supersymmetries of the action if they were to be interpreted as 
a "finite renormalization" of the superpotential. It then became clear (see [20] for a 
review and references therein) that the correct interpretation of these contributions was 
to consider them as IR singular D-terms in superspace, which are absent for instance in 
the more suitable Wilsonian definition of the effective superpotential. In what follows 
we would like to show that also in the case of ABJM theories does exist a special limit 
of vanishing external momenta for the vertex function which gives rise to a finite result. 

The vertex function, with the IR safe gauge choice, is guaranteed to be finite as 
long as the momenta are non-exceptional. By exceptional we mean when there exists 
at least one equation of the form ^ . piPi = with pi either or 1 and not all nor all 
1. In our case, it is easy to see that many exceptional configurations produce spurious 
IR divergences, for example if we choose any of the four momenta, say pi to be zero. 

If we were interested in finding an exceptional configuration which is IR safe and 
which leads to a constant, we would need at least two supplementary "exceptional" 
equations. In fact, we found that modulo equivalent choices, there is only one such 
choice of exceptional momenta which is IR finite. This is given by choosing Pi+ P2 = ^ 
and Pi + P4 = 0. We proceed to evaluate the graphs for this choice. 
For graph b we obtain 




2(P3+P4)' 




(3.3i; 



The reader might be worried that we put this graph to zero in the exceptional config- 
uration because of the (pa + ^4)^ = (pi + ^2)^ numerator without taking into account 
that the integral multiplying it is infrared divergent when P3+P4, = 0. A careful power 
expansion in |p3 +P4I gives 



2(P3+P4)' 




1 ^ 




Pi 


167r 







log(|) 



IGtt^ pI 



pi 



(3.32) 

where IC{z) is the complete elliptic integral of the first kind^ and the ellipsis are for 
higher orders in |p3 + ^4]. From this equation we see that ^3 + P4 — )■ is well defined 
and zero since all the coefficients in the expansion are finite in this limit. 



Notice that the coefficient in front of \p3 + Pi\ is imphcitly symmetric under o p4 due to 
71 r / 7:2-1 

In fact, all the coefficients of the expansion have this 



the property K. 
symmetry. 



1 



Pi 



Mir 

iP3n 



1-% 
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Graphs a may be represented in terms of elementary and Mellin-Barnes integral 



functions of the Lorentz invariants x 



(P3+P4)2 



and y 



Pi 



(P3+P4)2 



given by 



V 




V 





+ ^ d^' - ^^2(1 -x)- Li2{l -y)- log(x) log(y)) 

ioo 

|,.,,r-(-,)r(i-,)rMy^r,*^,)r(i^* + .),.. 

■ioo 

(3.33) 



+ 



647r2 
/? 1 



647r3 {2my 



This expression admits a well defined limit for (pa +^4)^ — given by 



327r2 



arccos 



b4 



arccos I - — r 1 + - log I -| 



(3.34) 



If we consider more in particular that ps + P4 = 0, then not only (ps + p^Y = but 
also P3 = P4, we find 



V 




V 




0. 



(3.35) 



Now we move on to graphs d. By expanding the products of momenta in the 
numerator of the integrals and properly completing squares (see trace properties of 7^^ 
matrices in the appendix), one can compare the resulting expression with the squared- 
completed expression of graphs a to conclude that 



V 




V 





V 



2(P3+P4)' 



V 





1 

32' 



Thus, according to the analysis we made before, in the limit pi + p2 
we obtain 

1 

~^ ~32' 



V 




+ v 




(3.36) 

-P3 - P4 ^ 

(3.37) 



^Definitions, properties and relevant references of Mellin-Barnes representation are given in the 
Appendix C. 
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Finally, it is possible to calculate the Feynman integral of graph e when pi +P2 = 
and Pi + P4 = by substituting pi = p, P2 = —p, Pa = —p and from momentum 
conservation p3 = p to obtain 



With all these elements we may make the sum to find the finite two-loop contribution 



T^^\A,B] = \\-8 - In') [ cfOd'x ^ e^ce'^'^tr [B^'AdB^Ab) • (3.39) 



As mentioned before we expect this to be a well defined and gauge invariant result. 
Nevertheless, it' s easy to show that if such contribution had to be interpreted as 
a "finite renormalization" of the superpotential it would inevitably break extended 
supersymmetry. 

4. General analysis 

We would like now to make some comments on the generality of our results. In a three 
dimensional theory there are two sources of infrared divergences in Feynman integrals. 
On the one hand we have the insertion of self energy corrected lines which may produce 
high powers of the propagators with a > |. On the other hand the presence of 
a three-lined vertex interaction with no external legs is potentially dangerous since, if 
there are only scalar propagators attached to it (no momenta in the numerator), an IR 
divergence is produced after loop integration. 

In general M = 2 Chern-Simons-Matter theories there are three-lined vertexes that 
couple chiral fields with the gauge vector and there is also the three gluon vertex. 
Consider the matter- vector coupling as shown in figure (6). 




2 




Figure 6: Matter-Gluon coupling 
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If wc integrate by parts on vertex 4 at least one of the D-operators of the gluon 
propagator, we get 

j2—^i^,^) -p-<^(4,2) ^-,2 ^(4,3) = ^2/2(fc_/)2 ^a()(l,4) D^D 5(4,2) D D 5(4,3)- 

(4.1) 

The appearance in the numerator of one of the momenta carried by the lines eliminates 
the IR threat as long as no self energy corrections are involved in the full graph (we 
deal with them in what follows). A similar analysis can be done for the three gluon 
vertex. 

It is quite obvious that the insertion of self-energy matter corrected lines inside 
any given graph, does not lower the scaling of the propagator thus not leading to IR 
issues. Then we conclude that IR problems are only generated by the insertion of self- 
energy corrected gluon lines: with the aid of the modified propagator we proposed in 
the introduction, it seems plausible that IR divergences can in principle be cured to all 
loop orders. 

To leading order we have shown that the key in the elimination of IR divergences 
was the completion of the 1-loop corrected gauge vector by adding the longitudinal part 
with the rj piece of the propagator. Having understood this mechanism that improves 
the IR behaviour of the gauge propagator correction at 1-loop, we may generalize this 
notion to all orders in A. Due to gauge invariance and parity we know [12] that the 
all order IPI vector self-energy calculated with the ordinary piece of the propagator is 
given by 

(4.2) 

where the coefficients Ai{k) and Bf{k) are functions that contain e-powers of the mo- 
mentum. That is, odd loop corrections contain the superspin 1/2 projector, and even 
loop corrections reproduce the original structure of the action. Any odd-loop correction 
from (4.2), when attached inside a graph produces a propagator given by 

-J2Bmx'^^^'^5,e,o,, (4.3) 

1=0 ' ' 

which, as noted before, will produce an IR divergence. On the other hand even-loop 
corrections, when attached inside a graph, produce a term ~ D°'Da/k'^ which behaves 
in the same way as the basic propagator, thus not leading to IR issues. These formulas 
can be readily derived using (A. 8). 



- 18 - 



Having understood the effect of corrected vector propagator insertions in graphs, 
we can now proceed to fix the rj parameter perturbatively as an odd power series in 
A. After fixing it to order one, 77^^^-) = — GttG^j^ ^^{p'^)~^X + O(A^), one calculates every 
connected (not only IPI) self energy vector correction at order A^, including lower loop 
O(A^) corrections with the ^77^^ piece of the propagator. With this result we fix the next 
coefficient ^77^^^ such that we complete the transverse projector with the longitudinal 
one effectively removing the source of infrared divergence at order A^. This process 
may be continued recursively thus improving the IR behavior of the propagator to all 
loops. 

In this way, if one considers a given graph which contains an L-loop-dressed gauge 
vector, then if L is odd there will always be a complementary graph in which we 
substitute that dressed line with the 77 piece of the propagator at the corresponding 
order in A, such that the whole line will behave as ~ instead, when L is even, 

the line behaves as the ordinary propagator ~ ^ -^^V.^') g^^^j needs no modifications. 
In both cases the graph will be IR safe. 

Some comments are in order about the ultraviolet behaviour of the non-locally 
gauged fixed theory. When we studied the effect of the 77 insertion in the case of matter 
self-energy graphs and superpotential corrections, we saw that the renormalization 
properties were not modified up to two loops. This makes this alternative gauge fixing 
procedure consistent, since we expect any gauge independent quantity of the theory 
to be independent from the procedure. As a further non trivial check, we have also 
verified that the renormalization properties of the theory in the gauge vector sector are 
also not modified up to two-loop order. 

5. Conclusions 

We studied the infrared behaviour of the off-shell amphtudes in three-dimensional 
Chern-Simons-matter theories with specific attention to the ABJM model. In J\f — 2 
superspace IR divergences show up in a very similar way as in four-dimensional Super- 
Yang-Mills theory, being related to the corrected vector superfield propagator inser- 
tions. At first, we showed that if the theory is gauge fixed in a standard fashion there 
is no way to get rid of the divergent integrals without losing the hermiticity of the ac- 
tion. Then we introduced a non-local gauge fixing procedure which leads to divergences 
cancelation without spoiling the renormalizability of the theory. In order to do so, the 
gauge-fixing parameter had to be perturbatively fine tuned. Moreover, we found in our 
computations that infrared infinities seem to be always associated to gauge dependent 
parts in the amplitudes, thus not affecting the physical quantities of the theory. As a 
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non-trivial output of our calculations we provided the two-loop finite correction to the 
IPI vertex function for ABJM theory in equation (3.39). 

It would be interesting to address the same problems in the case of Chern-Simons- 
matter theories described in Af — 1 superspace, starting for instance from the formula- 
tion of BLG theory [13] given in [14]. In this case the analogy with the four-dimensional 
case is lost and one might expect a different infrared behaviour of superspace propaga- 
tors. It would be also interesting to perform a similar analysis in the J\f — 3 harmonic 
superspace formulation of [18]. 
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A. Superspace Notations 



We use three-dimensional superspace notations adapted from [7] . We work in Euclidean 
space with a trivial metric 77^^ = 77^,^ = diag(l, 1, 1) and with Dirac Matrices {'y^)£ — 
iici, (72, (T3). We raise and lower spinor indexes through ip"' = C^^ipp and ipa = ip^C^on 
where the antisymmetric symbols C"^ and Cap arc defined by C^^ = — C12 = Notice 
that with this convention contractions arc always made going from the upper left corner 
to the lower right corner such that Grassmannian bilinears do not pick a sign after 
hermitian conjugation: (V"^)^ = (V'^V'a)^ = ^"V'a = Gamma matrices satisfy 

from which one can derive many useful trace properties such as 

rr(7'^7''7^7'') = 2r]^''r]P'' - 2r]''''r]'"' + 2r]^''rfP. (A.2) 

The similarity between H—1 d — A superspace and M—2 d=3 superspace is based on 
the isometry between both isospin groups: U{1) for (i = 4 and S0{2) for (i = 3. M=2 
superspace has two 2-component anticommuting coordinates 9'^ and 9^ such that if one 
defines complex coordinates 9" — 9^ — i92i 9°^ = 9"+i92 and complex spinor derivatives 

de. = ^(aw + da = ^(aw - (a.s) 

they satisfy da9^ — 5^, dofi^ — 5f , 8^9^ — 0, dofi^ — 0. Govariant derivatives are 

Da = \{D^^^ + iD^:'^) ^da + \9^dap, = \{D^^^ - iD^^^) = 9. + (A.4) 

such that they carry a representation of the super-algebra 

{Da,Dp}^i{-i^')apd^ = idap, = 0, = 0. (A.5) 

Apart from the fact that one does not make distinctions between dotted and un-dotted 
spinor indexes, this is the same algebra of covariant derivatives of A/'= 1 (i = 4 superspace 
thus making Feynman supergraph rules very similar to the known rules. On the other 
hand, one may construct contractions that were not allowed in four dimensions such 
as D°'Da or 9°'9a- An important property of the vector representation is that it is 
symmetric: C^^Pa^ — and C^^dap — 0, which is evident after one realizes that 7'* 
matrices with both spinor indexes up or down are symmetric with respect to those 
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indexes. Defining □ = d^^d^, = ^d'^^do.p, = iD'^Do, and = iD'^Do, the following 
properties hold 

D^D^ = 0, D^D^ = 0, [D", D^] = iO^^D^, [D^, D^] = iO'^^D^ 



jjijjijyL = UD^, D'^Df^ = D'^Dp = (A.6) 



Superspin projectors are defined as 

Vo = ^{D'D' + D'D'), V„2 = -^D^D^D^, (A.7) 
and together with D°'Da operator, they satisfy the useful properties 

{D'^D^f = nV,/2, Vi,2D''D^ = D^D^, V^D'^D^ = 0. (A.8) 

Our conventions for integration are J d?6 = | / dO^dOa, f d?0 = | / dO'^dOa and J d^O = 
J d'^9d'^9, such that up to a total space-time derivative 



d'e... = D\.. \e=-e=Q and J d'O . . . = D\ . . \e=e=o- (A.9) 
Finally, we define the 0-space 5-function as 5\e - 6') = {9- e'f{e - O'f. 



B. Feynman Rules 

We list some of the Feynman rules for ABJM theory. The vector superfield propagators 
are given in the a-gauge by: 

^"^^AA,^'''= J2 (^"^" + + "^') ^tofi'/c^t 
while using the ?7-gauge we obtain: 



ya yc f D'^D^ VJj)) \ 4 ^ . 



ya yc f D^D^ 7],{p) \ ^ 
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From the ghost action in (2.11) we find the ghost propagators: 



1 



K ....^...P'd = -±^S\9 - 9')S'X, ....y..fd = ^S\9 - d')5\5i (B.3) 



and from S^at we obtain the matter field propagators 

Ba ^ B^= ^^6\e - e')6i 51 si A^^A^= ^^5\9 - e')Si 6t 4, (B.4) 

where lowercase gauge indexes are omitted if they correspond to the same uppercase 
fiavor index {e.g. Ba = {Ba)^ -^^ = 

Some of the non trivial color structures of the vertexes we will need are 

B^\ /Ab 

X = ^^^Ace^"" i^l n Si - SI St Si) , (B.5) 



V 



c 




e^^esi. {si Sl 5S S^ - sf s\ 5^ St) , (B.6) 



Be ^ 



B^ = \^Si Si [st 5^ Si + S^, Si si) , (B.7) 



{Sl Si si - sl si 5f) [D'^{q)D^{p) - D-iq)D^ip)] , (B.8) 



and similarly for other vertexes involving A and V fields. We apply the usual D-algebra 
rules and regularize integrals when needed using dimensional reduction prescriptions. 

C. Relevant Integrals 

In the computation of two-point functions we introduced the function defined as 

f d'^k 1 _ G^a.fe) T{a + b- d/2)T{d/2- a)T{d/2-b) 



{27rY P« {k + pf^ ^p2y+b-d/2 ^4T,)d/2Y^a)r{b)r{d - a - b) {p^) 



2\a+b-d/2 



, (C.l) 
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and a particular two loop infrared divergent integral 

= I (2nf^p^i + ky{k + pfk^ = ^(1,3/2+.) (p')-'^ (C.2) 



All the relevant integrals in the calculation of the propagators can be reduced to the 
G-function form. 

In the computation of four-point integrals in the exceptional configuration, we 
found it necessary to expand Feynman integrals in powers of the kinematic invariants 
in order to carefully take the appropriate limits. In the end, the correct limit through 
this analysis became a confirmation of the "naive" result since it coincides with the 
limit taken directly on the integrand. To deal with the expansions we used multiple 
Mellin-Barnes contour representation of vertex integrals [15, 16]. These representations 
are based on the identity 

dsT{-s)T{s-ra)\^^ , (C.4) 



(A;2 + M2)« {M'^YT{a)2m J v / v ■ ' ^m"^ J ' 

— ioo 

where the contour is given by a straight line along the imaginary axis such that inden- 
tations are used if necessary in order to leave the series of poles s = 0, 1, ■ ■ ■ , n to the 
right of the contour and the series s = —a, —a — 1, ■ ■ ■ , —a — n to the left of the contour. 
After Feynman-parametrizing a triangle integral and using (C.4), the following formula 
holds 

r d<^k 1 (47r)-'^/^(p + g)^W^-^^^^) 

J {2'KYk'''^^{k-pYi^^{k + qyi^^~ ^^^(//^)^(rf-E^/^^) 

2 



X 



{p + q)y \{p + qy 

and for a vector-like triangle we have 

d'^k r (47r)-<^/2(p + g)2W2-EiMi^ / ^ ^2 



/ \27rY F/^i {k - pYi^- {k + qff^^ V{^Xi)V{d - Ei /^^ + 1) V27ri 

ioo 

X j dsdtT{-s)V{-t)T{i,, + s + mY,p,-i + s + t)(^j-^^ 



X 



2 \ * / 2 \ i 



X 
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[r(f-/xi-/x2-s)r(f-/xi-/x3-t+i)p"-r(f-/xi-/x2-s+i)r(f-/xi-/x3-t)g"] , 

(C.6) 

where the multiple contours are taken using the same convention as the first definition 
unless otherwise indicated. It is customary to indicate with a * over the r{z) function 
the case where one leaves a pole to the other of the conventional side of the contour. 
With these representations, among with Barnes 1** and 2"'' lemmas 

1 f J ^T^ru NT./ NT./ J N r{a + c)r{a + d)r{b + c)r{b + d) 
1 f r(a + s)r(b + g)r(c + s)r{d - s)r(-s) 

with e = a + b + c + d, (C.8) 



2ni J r(e + s) 

r{a)T{b)T{c)T{a + d)T{h + d)T{c + d) 
r(e - a)r(e - b)T{e - c) ' 



and their multiple corollaries, we were able to carefully expand the 2-loop four-point 
integrals in the relevant kinematic invariants in order to take the limit of exceptional 
momenta. 
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